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ABSTRACT 


Since the Voyager 1 (VI) flyby in 1980, Titan’s exploration from space and the ground has been on- 
going for more than a full revolution of Saturn around the Sun (one Titan year or 29.5 Earth years was 
completed in May 2010). In this study we search for temporal variations affecting Titan’s atmospheric 
thermal and chemical structure within that year. We process Cassini/CIRS data taken during the Titan 
flybys from 2006-2013 and compare them to the 1980 Vl/IRIS spectra (re-analyzed here). We also 
consider data from Earth-based and -orbiting observatories (such as from the ISO, re-visited). When we 
compare the CIRS 2010 and the IRIS data we find limited inter-annual variations, below the 25% or 
35% levels for the lower and middle, or the high latitudes, respectively. A return to the 1980 
stratospheric temperatures and abundances is generally achieved from 50°N to 50°S, indicative of the 
solar radiation being the dominating energy source at 10 AU, as for the Earth, as predicted by GCM 
and photochemical models. However, some exceptions exist among the most complex hydrocarbons 
(C 4 H 2 and C 3 H 4 ), especially in the North. In the Southern latitudes, since 2012, we see a trend for an 
increase of several trace gases, possibly indicative of a seasonal atmospheric reversal. At the Northern 
latitudes we found enhanced abundances around the period of the northern spring equinox in mid-2009 
(as in Bampasidis et al. 2012), which subsequently decreased (from 2010-2012) returning to values 
similar to those found in the V 1 epoch a Titanian year before. 

Key words: infrared: solar system - planets and satellites: individual (Titan) - atmospheres - planets 
and satellites: composition - radiation mechanisms: thermal - radiative transfer 


1. INTRODUCTION 

Mid-2010 Cassini observations marked one Titan year since the first extended space exploration 
of the satellite by Voyager 1 (VI) during the 12 November 1980 flyby. Titan, Saturn’s largest moon, 
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follows its planet on its course around the Sun on an inclined orbit (26.7°) which gives rise to seasons, 
similarly to Earth. The Cassini mission had performed 68 Titan flybys by May 2010 (See Tables 1-4 in 
Bampasidis et al. 2012) when the 9° solar longitude (Ls), was almost exactly the same as at the time of 
the VI flyby (just after spring equinox in the northern hemisphere). Since then, and until the end of 
2012, more mid-infrared nadir spectra of Titan were acquired, which we analyze here. These data allow 
us to probe Saturn’s and, especially in what concerns us here, Titan’s stratosphere (i.e. altitudes 
between -100 and 500 km above the surface), where neutral species include hydrocarbons, nitriles and 
oxygen compounds. 

Stratospheric temperatures and the trace gaseous content, as well as the associated vertical and 
meridional (latitudinal) distributions, have been reported both from Voyager and Cassini data 
(Coustenis and Bezard 1995; Achterberg et al. 2008, 201 1; Coustenis et al. 1991, 2007, 2010; Teanby 
et al. 2008, 2009, 2012; Vinatier et al. 2007, 2010a, b). Here we search for temporal variations within 
these structures in Titan’s stratosphere. Such an investigation was first performed between the two 
Titan Voyager 1 and 2 thermal IR datasets (Letourneur and Coustenis 1993) just months apart (the 
results were quite similar within error bars). Since then, other investigators have reported seasonal 
variations and trends during the Cassini mission (e. g. Teanby et al. 2008, 2009, 2010, 2012; Jennings 
et al. 2012; Bampasidis et al. 2012). 


2. OBSERVATIONS 

The Vl/IRIS instrument, a Michelson interferometer that covered the 180-1500 cm' 1 range 
(Hanel et al. 1980; Flasar et al. 2004), acquired about 300 usable spectra with a resolution of 4.3 cm' 1 
during the 1980 Titan flyby. Its heir, the Composite Infrared Spectrometer (CIRS), acquired a large 
number of spectra during the Cassini Titan flybys (about 90 by the end of 2012), in the thermal infrared 
range from 600 to 1430 cm' 1 (focal plane 3, or FP3 (600-1 100 cm' 1 ), and 4 or FP4 (1100-1430 cm' 1 ) in 
nadir and limb geometry conditions at high, medium and low spectral resolutions (0.53, 2.54 and 15 
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cm’ 1 respectively; Flasar et al. 2005). The CIRS data represent a significant improvement in number of 
spectra, spatial coverage, data quality and spectral resolution (by about an order of magnitude, see 
Figure 1). 

We analyze several averages of CIRS nadir-viewing FP3 and FP4 spectra at medium and high 
resolution during the Cassini mission from July 2004 until 2013, complementing our previous analysis, 
which stopped in 2010 at flyby T72 (Bampasidis et ah, 2012), by including data from 2011 and 2012 
Cassini/CIRS spectra at 50°N, 50°S and near-equator. In particular, we have added one selection for 
50°N recorded in September 2012 (T86 flyby) with 761 spectra (35.5 S/N and 1.24 airmass) in FP3 and 
831 spectra (57.6 S/N and 1.28 airmass) in the associated FP4. Similarly, for 50°S a selection from 
September 2012 contains 777 spectra (22.5 S/N and 1.34 airmass) in FP3 and 612 spectra (30.4 S/N 
and 1.27 airmass) in FP4. At the equatorial latitudes, we have included three more selections recorded 
in January, February and May 2012 during the T81, T82 and T83 Cassini flybys respectively, with 
1431, 1147 and 617 spectra for each FP3 (S/N : 57.7, 53.3 and 34.5 ; airmass : 1 . 1 9, 1.01 and 1.01, 
respectively). In the FP4 region we have 1469, 1131 and 551 spectra for each flyby with S/N of 129.5, 
116.5 and 79.9 and 1.10, 1.11 and 1.01 airmasses, respectively. 

We have also compiled a listing of stratospheric-pertaining ground-based and Infrared Space 
Observatory (ISO, Coustenis et al. 2003) abundance inferences. These values are not disk-resolved and 
some times refer to different levels in the atmosphere of Titan, so that direct comparisons are difficult 
with the space data, all the more so since the error bars are not firmly defined in most of the ground- 
based studies. These refer mostly to molecules such as C 2 H 6 and HCN (Table 1, machine readable) and 
the snapshots we have between 1980 and the Cassini mission show no large departures from one season 
to the other for these few molecules. The data from ISO/SWS taken in 1997 concern more molecules 
and are briefly discussed hereafter, where we mainly focus on comparisons between CIRS and IRIS, as 
these space datasets are the most directly comparable. 
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3 . ANALYSIS 


In this work we use Cassini/CIRS data to extend the 2006-2010 study of Bampasidis et al. (2012) 
until the end of 2012. We combine this with a re-analysis of the Voyager data and ISO datasets in order 
to define the neutral chemistry and atmospheric evolution on Titan on timescales longer than a Titan 
year. We focus on latitudes from 50°S to 50°N as these have the most complete temporal coverage. 

We use our radiative-transfer code ARTT (Coustenis et al. 2010; Bampasidis et al. 2012), which 
includes several new molecules and the updated spectroscopic parameters for the gases and their 
isotopologues (Coustenis et al. 2003, 2008; Jennings et al. 2008; Nixon et al. 2008a, b, 2009; Rothman 
et al. 2009; Jacquinet-Husson et al. 201 1). The thennal structure is derived using an a priori approach 
described in Achterberg et al. (2008, 2011) and applied on the Cassini/CIRS and Vl/IRIS data by 
fitting the V4 methane band at 1305 cm' 1 (Figure 1). The altitude/pressure ranges probed by the 
temperature profiles depend on the latitude, but are in general between 0.1 and 10 mbar (for details see 
Coustenis et al. 1995, 2003, 2010; Bampasidis et al., 2012). We then use the temperature profile in our 
radiative transfer code to derive the chemical composition of the trace gases in the rest of the spectrum. 

We apply this method to different latitudes from 50°N to 50°N for CIRS and IRIS. The mixing 
ratios are assumed to be vertically-constant profiles above the condensation level and the altitudes 
probed for all space data correspond in general essentially to 100-250 km (0.5-10 mbar, see Coustenis 
et al. 2010 and Bampasidis et al. 2012 for more details). When the C2H2 band is not properly fitted at 
high or low resolution in the entire range with such a profile, we give an averaged abundance over the 
700-750 cm' 1 range in the results (for details see Coustenis et al. 1995, 2010). The contribution 
functions for these molecules remain within the ranges detennined in these previous papers. 

Our comparison with Vl/IRIS date is restrained to those molecules with strong enough emission 
bands to be unambiguously detected also in the Voyager spectra with high enough S/N ratio and 
sufficient accuracy. These molecules are: CO2, HCN, C2H2, C2H6, C3H4, C4H2. In this analysis we do 
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not include the results we infer for C3H8 and C2H4, which suffer from higher uncertainty than the other 
molecules, since their weak and/or blended bands are difficult to exploit in the Vl/IRIS spectra. 

Table 2 shows the re-analyses of the Vl/IRIS data (updating Coustenis and Bezard 1995) and of 
the ISO Short Wavelength Spectrometer (SWS) disk-averaged spectrum (updating Coustenis et al. 
2003) using the same model as for the CIRS data. The differences with the previous publications are 
minor, except for the molecules that benefited from significantly improved spectroscopic parameters. 

4. RESULTS ON THE STRATOSPHERIC STRUCTURE DURING ONE TITANIAN 

YEAR 

The best fits obtained for the Vl/IRIS and the January 2010 Cassini/CIRS 50°N selections are 
shown in Figure 1. 

4. 1 Inter-annual temperature variations 

We compare the new temperature profiles inferred for the Voyager selections with the ones found 
for 2010 from CIRS data at 50°N, equator and 50°S (Fig. 2). We also overplot on that figure the 
temperature profile that corresponds to the disk-average ISO observations of 1997, which pertain morst 
closely to conditions at mid-latitudes. 

I 11 the southern hemisphere, the 1980 and 2010 temperature profiles are quite similar, with very 
small differences of less than 5 K and more typically 2-3 K between 0.1 (250 km) and 10 mbar (100 
km), in the region probed by our data. Above the 0.1-mbar pressure level the temperature appears to be 
cooling, e.g., by as much as 15 K at 0.01 mbar (400 km) from 2005 to 2010, but these inferences are 
very uncertain because the nadir data do not probe such high altitudes. The reader is referred to limb- 
data analyses as, for instance, in Vinatier et al. (2007, 2010) for the temperatures at higher altitudes. 
The revised Voyager temperature profiles are then consistent with the CIRS 2010 ones in the lower 
latitudes and no significant variations are detected. 
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Similarly, no significant change was found in the equatorial mid-stratospheric temperatures, with 
the 1980 and the 2010 temperature profiles always within 5-6 K of each other up to the 0.1-mbar level, 
although again it appears that a cooler equatorial mesosphere is observed after one Titan year as in the 
South. We note that the temperature profile corresponding to the ISO 1997 data follows very closely 
the CIRS 2010 profile throughout the equatorial stratosphere. All in all, the average vertical 
stratospheric temperature profiles from 50°S to the equator appear quite to each other at all times. The 
50°S temperatures are about 3-6 K cooler than the equatorial ones for CIRS (see e. g. Bampasidis et al. 
2012 for more precise information on the thermal evolution during the Cassini mission). 

At northern latitudes, temporal variations were reported during the Cassini mission, while the 
shape of the thennal structure changed towards a vertically more homogeneous profile in the past 4 
years not showing the marked stratopause observed earlier (Achterberg et al. 2011, Bampasidis et al. 
2012). However, when we compare the 1980 and the 2010 stratospheric profiles at 50°N, we find the 
VI profile to be quite compatible to the CIRS one, equivalent within 3-4 K up to the 0.2-mbar (220 
km) level. Thus, in the stratosphere, the thennal structure has returned to the VI shape after a Titanian 
year (Fig. 2). 

4.2 Chemical composition variations 

With the temperature profiles determined as indicated above, we infer the chemical composition 
for Titan’s stratosphere from CIRS, IRIS and ISO spectra. The new and more accurate chemical 
composition inferred at the time of the Voyager 1 encounter (12 Nov. 1980) is given in Table 2 for the 
latitudes of interest here. The revised abundances for the ISO/SWS data of 1997 are also listed in that 
Table. 

The ratios of the CIRS calculations (of high resolution, 0.5 cm' 1 ) with respect to the Vl/IRIS 
values (dividing each CIRS result with the Vl/IRIS 1980 inference) are shown in Fig. 3 for 50°S, 
equator, and 50°N. Since the model parameters for all these calculations are the same, we only consider 
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the relative uncertainties that are generally small for the CIRS data (since the noise is small, given the 
large number of spectra summed) but mainly concern uncertainties as to some of the spectroscopic data 
and temperature profile inferences (see Bampasidis et al. 2012 for details). In Fig. 3 we display only 
these relative error bars for all the molecules so that we can derive the variations with time within the 
Cassini mission period for each species. Note that the ratios are not indicative of any variations wrt to 
V 1/IRIS except for the 2010 period, when they are meaningful. Indeed, for the time when we return to 
the VI epoch, one needs to consider in addition, when comparing CIRS and IRIS, the uncertainty due 
to the Vl/IRIS random errors (noise, calibration, temperature profile determination and lit of the 
continuum). When this is taken into account, at the 2010 dates indicated by arrows in Fig. 3, the level 
of uncertainty on the ratio (total between the IRIS and the CIRS error bars) is 27%, 23% and 35% for 
50°S, equator and 50°N, respectively. 

Hereafter we compare for each molecule: a) the trend during the Cassini mission shown by the 
ratios; b) the comparison with the abundances a Titan year ago, around mid-2010. 

4.2.1 Seasonal evolution of the stratospheric composition 

Starting with the lower latitudes (50°S), we find that during the Cassini mission, the CIRS high- 
resolution data indicate that, within measurement uncertainties, CCF, C 2 H 2 and C 2 H 6 can be considered 
unvarying to within 25% up to 2013, as indicated also in Bampasidis et al (2012), while HCN, C 4 H 2 
and C 3 H 4 show some excursions, especially in the more recent years (Fig. 3a). We further note that the 
results from the analysis of the Sept. 2012 selection indicate that there is a significant increase of the 
abundances of some of the molecules like C 4 H 2 , C 3 H 4 , CO 2 , and HCN, while other more abundant 
molecules (C 2 H 2 and C 2 H 6 ) remain constant. This increase in abundance is compatible with the 
expected seasonal reversal in Titan’s stratospheric composition for some gases and the haze (see 
Teanby et al. 2012; Jennings et al., 2012 and references therein). 
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Fig. 3b shows the situation for the equator during the Cassini mission and until 2013, in which 
we find, as for 50°S, rather constant variations with time for all molecules. C 4 H 2 appears to have 
decreased in abundance from 2008 to 2012 by 30%, but the abundance of this molecule suffers of large 
uncertainties due to its lacking spectroscopic parameters. Since 2011, the concentrations of all the gases 
show a trend for a certain decrease, which will need further confirmation. 

As noted in Bampasidis et al (2012), at 50°N, we also find that, with the possible exception of 
C 2 H 6 , CO 2 and C 3 H 4 , all the molecules show significant variations with time during the Cassini 
mission. Thus, while from 2007 to mid-2008 the gas mixing ratios are constant to within 30%, we 
detect a general trend for a pronounced increase from mid-2008 to mid-2009 with a maximum (by 
factors of 2-3) for almost all molecules in April-May 2009 (the closest to Titan’s northern spring 
equinox, NSE, August 2009). This is followed by a strong decrease, which reduces the 2010 
concentrations to the values of the period preceding the increase (Fig. 3c), as confirmed by the Sept. 
2012 data. 

4.2.2 Inter-annual compositional variations 

We compare CIRS 2010 abundances with the retrievals of the V 1/IRIS (1980) re-analyses by 
computing the ratios for the observed molecules. For comparison’s sake with the Vl/IRIS data, we 
have included some abundance retrievals from medium-resolution (2.5 cm' 1 ) CIRS spectra for 2010. 
The arrows in Fig. 3 point to the V 1 encounter. 

At 50°S, the 2010 April high-resolution and May- June medium-resolution spectra are very 
compatible with each other. However, the abundances derived from the medium-resolution spectra are 
generally slightly higher by 10-15% (Fig. 3a). All gaseous abundances from the CIRS datasets are 
within ±25-30% of the values observed in 1980 by Voyager 1 for both the higher and lower resolution 
measurements, which, given the difference in resolution and attitudes probed, can be considered as 
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quite compatible. Exceptions are the more complex hydrocarbons (C 3 H 4 and C 4 H 2 ), which are below 
the Vl/IRIS abundances at about 30-40% (Fig. 3a). 

The 2010 equatorial CIRS mixing ratios of all the species from all resolution data are within the 
23% uncertainty, except for C 4 H 2 , which is again about 35% lower than Vl/IRIS. 

In general, at 50°N we find in 2010 a return to the Vl/IRIS values within the 35% uncertainty 
level for all molecules except the two complex hydrocarbons C4H2 and C3H4, which stay well below 
the 1980 abundances by about 50-55%. Again, the May 2010 50°N CIRS nadir data at 2.5 cm " 1 are 
close to the high resolution values (differing by no more than ±15%) as also noted in Bampasidis et al. 
(2012) (compare the open and full circles on Fig. 3c for 2010 values). 

We also considered ground-based and disk-averaged ISO data obtained before Cassini or more 
recently with Herschel (see Table 1). These data are limited to a few molecules and do not show, as far 
as we can see, any significant variations within a Titan year. Only the ISO data, taken toward the 
northern fall equinox (NFE, occurring in 1995), yield somewhat higher values for the species discussed 
here, such as HCN. One of the rare models dealing with long-tenn variations in Titan’s atmosphere 
over a Titan year starting in 2000 (Hourdin et al. 2004, their Fig. 9), shows that soon after NFE, in the 
80-240 km altitude range, the high-latitude HCN concentration grows very fast under the effect of 
downward advection, as do the latitudinal contrasts with low latitudes, reaching a maximum, so that the 
ISO disk-average data, which include some of the polar enhancement in gaseous components observed 
at NFE, may indeed exhibit higher values than any disk-resolved space equivalent data. However, it is 
very difficult to conclude that the ISO values are indeed higher than the Vl/IRIS or the CIRS ones, 
given the impossibility to properly assess and compare ground-based or Earth-bound observations with 
space in situ ones given the differences in altitudes and regions probed, in spectral resolution, in 
modeling parameters, etc. We note, however, that the HCN stratospheric vertical profile obtained by 
Herschel/SPIRE recently (Courtin et al. 2011) for July 2010 is in good agreement with our findings at 
around 5 mbar (125 km in altitude) for equatorial latitudes (~1.5 10' 7 , Fig. 3 in Bampasidis et al. 2012; 
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Fig. 5b and Table 2 of Coustenis et al. 2010), given that we only consider constant vertical profiles. We 
deplore, in general, the lack of more ground-based observations of Titan’s atmospheric chemistry, 
which could provide good comparison points even on a disk-average level during a Titan year. 

We also note that the results of a recent analysis of CIRS limb data (Vinatier et al., 2010) at 50°S, 
1°N and 51°N data taken in June 2010, but covering higher altitudes (up to 500 km) and the simple 
vertical distributions of some hydrocarbons and nitriles retrieved from scarce horizontal viewing 
Vl/IRIS data (Coustenis et al. 1991) are compatible in the general trend showing the nitriles increasing 
with altitude faster than the hydrocarbons. But again, the comparison cannot be pushed too far given 
the differences in vertical resolution, modeling, gaseous distributions, thermal structure, etc. 

5. DISCUSSION 

The fulfillment of one Titan year of observations provides us for the first time with the 
opportunity to evaluate the relative role of different physical processes in the long-term evolution of 
this complex environment. By comparing VI (1980), ISO (1997) and Cassini (2010) data we find that 
the temperature structure is quite similar within error bars (Fig. 2). The chemistry also appears to be 
generally unchanging (Vl/IRIS vs. Cassini/CIRS are within 25-35%), with some exceptions, among 
the more complex hydrocarbons such as C3H4 and C4H2. 

The return of today's abundances to values close to the ones derived from Voyager spectra (at the 
same season) is an indication that, as for the Earth, the solar radiation dominates over other energy 
sources, even at 10AU. Nevertheless, the differences observed (lower C3H4 and C4H2 mixing ratios) 
and the slightly higher values observed from the ground and from ISO near the NFE could suggest that 
other processes may be at play as well, for example, the variability of the solar insolation itself through 
the 1 1-year solar cycle or complex circulation phenomena. We also note from the variations during the 
Cassini mission found for the gases investigated here, that in the South pole, some of the minor species 
tend to show a strong increase in 2012, in accordance with what is expected from seasonal models of 
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Titan, indicating that a reversal of the north-south asymmetry is on-setting, and confirming the trend 
seen at higher altitudes in CIRS limb data (Teanby et al. 2012). 

The detection of long-tenn temporal variations of the atmospheric chemical composition and 
thermal structure informs us on the impact of different processes such as photochemistry and dynamics. 
The main energy source defining Titan's chemical composition is solar radiation, while energetic 
particles from Saturn's magnetosphere, as well as galactic cosmic rays have also a significant 
contribution. All these energy sources vary temporally in magnitude during the course of one Titan 
year. C 3 H 4 and C 4 H 2 are short-lived and low-abundance complex hydrocarbons in Titan’s stratosphere 
and hence more prone to inter-annual variations. The temporal and spatial (due to Titan's inclination) 
variations in the energy input to Titan's atmosphere is a driver for changes in the advection patterns, 
which in turn provide a stronger variability in the latitudinal abundances of photochemical species. 
Furthermore, both the changes in the energy field and on the dynamics affect the production and 
evolution of aerosols in Titan's atmosphere. The abundance peak corresponds to a period during which 
the observed collapse of the detached aerosol layer (West et al. 2011) suggests that the dynamics went 
through a rapid transition which should influence the gas distribution and can be tied to the atmospheric 
circulation reversal predicted by GCM numerical models (Hourdin et al. 2004; Rannou et al. 2005; 
Crespin et al. 2008; Lebonnois et al. 2012). At the same time, the variable north/south asymmetry in 
Titan's albedo and haze content reveals the consequences of seasonal changes in the aerosol properties 
deeper in the atmosphere (Jennings et al. 2012). 

Circulation and photochemical models must satisfy the constraints set by these new results, and 
further observational constraints will come from data acquired during the Cassini mission and up until 
the next summer solstice (2017). We note that we are still lacking disk-resolved data of the atmospheric 
composition for seasons from summer to winter solstice in the northern hemisphere. 
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Molecule 

50°S 

7°N 

50°N 

ISO 

disk-average 

c 2 h 2 

2.25±0.45 x 10 ' 6 

2.62±0.4 x 10 ' 6 

3.35±1.0 x 10 ' 6 

4.95±0.5 x 10 ' 6 

c 2 h 4 

1.14±0.45 x 10 ' 7 

1.13±0.3 x 10 ' 7 

2.55 ±1.22 x 10 ' 7 

9.33±0.4 x 10 ' 8 

c 2 h 6 

9.50±1.85 x 10 ' 6 

1.01±0.15 x 10 ' 5 

1.35 ±0.34 x 10 ' 5 

1.32±0.5 x 10 ' 5 

C3H4 

5.60±1.75 x 10 ' 9 

7.80±1.6 x 10 ' 9 

3.10 ±0.93 x 10 ' 8 

1.44±0.4 x 10 ' 8 

c 3 h 8 

1.40±0.6 x 10 ' 6 

1.60±2.2 x 10' 6 

2 . 1 0 ± 1 .0 x 10' 6 

4.76±0.8 x 10 ' 7 

c 4 h 2 
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HCN 
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2.05±0.2 x 10 ' 8 


Table 2. New values of Vl/IRIS abundance retrievals at 50°S, equator and 50°N with 3-sigma error 
bars including both systematic and random uncertainties. The ISO values were also re-calculated with a 
temperature profile matching the emission observed in the CH 4 V 4 band at 1304 cm ' 1 and the new 
spectroscopic data available for some molecules. The corresponding altitudes are the same as in 
Coustenis et al. (2003). The re-analysis of the Vl/IRIS data allowed us to detect in them for the first 
time HC 3 N and C 2 HD (see fit in Fig. 1), albeit given the low spectral resolution of the Voyager 

spectra, we were unable to derive accurate mixing ratios for these molecules and therefore they are 
excluded from this paper. The same pertains to C 3 H 8 and C 2 H 4 which similarly bear large uncertainties. 
The CIRS abundances from 2006 to 2013 are shown with respect to V 1/IRIS in Fig. 3. 


FIGURE CAPTIONS 


Figure 1. Upper panel: fit of the Vl/IRIS Nov. 1980 data of Titan taken at 50°N, with the updated code 
and spectroscopic parameters. The spectral resolution is 4.8 cm' 1 . The improved fit in the region 670- 
700 cm' 1 , with respect to the Coustenis and Bezard (1995) publication, is due to the inclusion of the 
C 2 HD and C 6 H 6 bands which were not identified at that time. Lower panel: fit of the 0.5 cm' 1 resolution 
CIRS spectrum taken in Jan. 2010. Left: FP3 region; right : FP4 region. 


Figure 2. Retrieved thennal profiles from CIRS and Vl/IRIS data (in full and dashed lines 
respectively) at 50°S, 50°N and equator for the 2010 nadir data and IRIS 1980 data. The CIRS 
temperature profiles correspond to April, July and January 2010 at 50°S, equator and 51°N 
respectively. The temperature profile corresponding to the ISO disk-averaged SWS data of 1997 is also 
shown in dotted black lines. The pressure/altitude ranges probed by these profiles vary with latitude, 
but are generally in the 0. 1-10 mbar region (see text). The uncertainty is about 1 K at 1 mbar and 5 K at 
5 mbar for the CIRS profiles and about 1 K more at each level for the IRIS profiles due to the higher 
noise. For ISO, the error bars are about 2 K at pressures lower than 1 mbar and at most 5 K at higher 
ones (altitudes lower than 200 km). 


Figure 3 a) variations of the abundances of gases in Titan’s stratosphere with respect to the 
Vl/IRIS values at 50 °S (normalized as 1) for: C 2 H 2 (black), C 2 H 6 (red), CO2 (brown), HCN (blue), 
C3H4 (green), C4H2 (magenta). The arrows indicate the date of the 1980 VI encounter. Larger open 
circles correspond to medium-resolution (2.5 cm' 1 ) data that we have in 2010. b) : same as a) but for 
the equator, c) : same as a) but for 50 °N. The vertical error bars are only from CIRS relative 
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uncertainties (see text). In mid-2010, at the time of the return to the Voyager encounter epoch 
(indicated by arrows), all gases show values within the uncertainty levels of 27%, 23% and 35% for 
50°S, equator and 50°N, respectively (not shown in the figure to avoid confusion), which take into 
account the Vl/IRIS errors also. Exceptions are C3H4 and C4H2 in some cases (see text). 
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